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There have been several popular reports of various groups exploiting the deterministic nature of the 
game of roulette for profit. Moreover, through its history the inherent determinism in the game of 
roulette has attracted the attention of many luminaries of chaos theory. In this paper we provide 
a short review of that history and then set out to determine to what extent that determinism can 
really be exploited for profit. To do this, we provide a very simple model for the motion of a 
roulette wheel and ball and demonstrate that knowledge of initial position, velocity and acceleration 
is sufficient to predict the outcome with adequate certainty to achieve a positive expected return. 
We describe two physically realisable systems to obtain this knowledge both incognito and in situ. 
The first system relies only on a mechanical count of rotation of the ball and the wheel to measure 
the relevant parameters. By applying this techniques to a standard casino-grade European roulette 
wheel we demonstrate an expected return of at least 18%, well above the —2.7% expected of a 
random bet. With a more sophisticated, albeit more intrusive, system (mounting a digital camera 
above the wheel) we demonstrate a range of systematic and statistically significant biases which can 
be exploited to provide an improved guess of the outcome. Finally, our analysis demonstrates that 
even a very slight slant in the roulette table leads to a very pronounced bias which could be further 
exploited to substantially enhance returns. 
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"No one can possibly win at roulette unless he 
steals money from the table when the croupier 
isn't looking" (Attributed to Albert Einstein in'^ 

Among the various gaming systems, both cur- 
rent and historical, roulette is uniquely de- 
terministic. Relatively simple laws of motion 
allow one, in principle, to forecast the path 
of the ball on the roulette wheel and to its 
final destination. Perhaps because of this ap- 
pealing deterministic nature, many notable 
figures from the early development of chaos 
theory have leant their hand to exploiting this 
determinism and undermining the presumed 
randomness of the outcome. In this paper 
we aim only to establish whether the deter- 
minism in this system really can be profitably 
exploited. We find that this is definitely pos- 
sible and propose several systems which could 
be used to gain an edge over the house in a 
game of roulette. While none of these sys- 
tems are optimal, they all demonstrate posi- 
tive expected return. 



^'Electronic mail: michael.small@uwa.edu.au 



I. A HISTORY OF ROULETTE 

The game of roulette has a long, glamorous, in- 
glorious history, and has been connected with sev- 
eral notable men of science. The origin of the game 
has been attributecP, perhaps erroneousljl^, to the 
mathematician Blaise PascaP. Despite the roulette 
wheel becoming a staple of probability theory, the 
alleged motivation for Pascal's interest in the de- 
vice was not solely to torment undergraduate stu- 
dents, but rather as part of a vain search for perpet- 
ual motion. Alternative stories have attributed the 
origin of the game to the ancient Chinese, a French 
monk or an Italian mathematiciarP'^l jn any case, 
the device was introduced to Parisian gamblers in 
the mid-eighteenth century to provide a fairer game 
than those currently in circulation. By the turn of 
the century, the game was popular and wide-spread. 
Its popularity bolstered by its apparent randomness 
and inherent (perceived) honesty. 

The game of roulette consists of a heavy wheel, 
machined and balanced to have very low friction and 
designed to spin for a relatively long time with a 



slowly decaying angular velocity. The wheel is spun 
in one direction, while a small ball is spun in the 
opposite direction on the rim of a fixed circularly in- 
clined surface surrounding and abutting the wheel. 
As the ball loses momentum it drops toward the 
wheel and eventually will come to rest in one of 37 
numbered pockets arranged around the outer edge of 
the spinning wheel. Various wagers can be made on 
which pocket, or group of pockets, the ball will even- 
tually fall into. It is accepted practise that, on a suc- 
cessful wager on a single pocket, the casino will pay 
35 to 1. Thus the expected return from a single wa- 
ger on a fair wheel is {35 + l)x^ + (-1) w -2.79P. 
In the long-run, the house will, naturally, win. In the 
eighteenth century the game was fair and consisted 
of only 36 pockets. Conversely, an American roulette 
wheel is even less fair and consists of 38 pockets. We 
consider the European, 37 pocket, version as this is 
of more immediate interest to us^. Figure [l] iUus- 
trates the general structure, as well as the layout of 
pockets, on a standard European roulette wheel. 

Despite many proposed "systems" there are only 
two profitable ways to play roulett^. One can ei- 
ther exploit an unbalanced wheel, or one can ex- 
ploit the inherently deterministic nature of the spin 
of both ball and wheel. Casinos will do their utmost 
to avoid the first type of exploit. The second exploit 
is possible because placing wagers on the outcome 
is traditionally permitted until some time after the 
ball and wheel are in motion. That is, one has an 
opportunity to observe the motion of both the ball 
and the wheel before placing a wager. 

The archetypal tale of the first type of exploit 
is that of a man by the name of Jagger (various 
sources refer to him as either William Jaggers or 
Joseph Jagger, or some permutation of these). Jag- 
ger, an English mechanic and amateur mathemati- 
cian, observed that slight mechanical imperfection 
in a roulette wheel could afford sufficient edge to 
provide for profitable play. According to one incar- 
nation of the tale, in 1873 he embarked for the casino 
of Monte Carlo with six hired assistants. Once there, 
he carefully logged the outcome of each spin of each 
of six roulette tables over a period of five weekJ^. 
Analysis of the data revealed that for each wheel 
there was a unique but systematic bias. Exploiting 
these weaknesses he gambled profitably for a week 
before the casino management shuffled the wheels 
between tables. This bought his winning streak to a 
sudden halt. However, he soon noted various distin- 
guishing features of the individual wheels and was 
able to follow them between tables, again winning 
consistently. Eventually the casino resorted to re- 
distributing the individual partitions between pock- 



ets. A popular account, pubhshed in 1925, claims 
he eventually came away with winnings of £65, OOCP. 
The success of this endeavour is one possible inspira- 
tion for the musical hall song "The Man Who Broke 
the Bank at Monte Carlo" although this is strongly 
disputecP. 

Similar feats have been repeated elsewhere. The 
noted statistician Karl Pearson provided a statistical 
analysis of roulette data, and found it to exhibit sub- 
stantial systematic bias. However, it appears that 
his analysis was based on flawed data from unscrupu- 
lous scribeJSl (apparently he had hired rather lazy 
journalists to collect the data). 

In 1947 irregularities were found, and exploited, 
by two students, Albert Hibbs and Roy Walford, 
from Chicago Universitjff^'^. Following this line 
of attack, S.N. Ethier provides a statistical frame- 
work by which one can test for irregularities in the 
observed outcome of a roulette wheeP^. A similar 
weakness had also been reported in Time magazine 
in 1951. In this case, the report described various 
syndicates of gamblers exploiting determinism in the 
roulette wheel in the Argentinean casino Mar del 
Plata during 194^^. The participants were colour- 
fully described as a Nazi sailor and various "fruit 
hucksters, waiters and farmers'^. 

The second type of exploit is more physical (that 
is, deterministic) than purely statistical and has con- 
sequently attracted the attention of several math- 
ematicians, physicists and engineers. One of the 
firsiPl was Henri Poincar^ in his seminal work Sci- 
ence and MethoM While ruminating on the nature 
of chance, and that a small change in initial condi- 
tion can lead to a large change in effect, Poincare il- 
lustrated his thinking with the example of a roulette 
wheel (albeit a slightly different design from the 
modern version) . He observed that a tiny change in 
initial velocity would change the final resting place 
of the wheel (in his model there was no ball) such 
that the wager on an either black or red (as in a 
modern wheel, the black and red pockets alternate) 
would correspondingly win or lose. He concluded 
by arguing that this determinism was not important 
in the game of roulette as the variation in initial 
force was tiny, and for any continuous distribution 
of initial velocities, the result would be the same: 
effectively random, with equal probability. He was 
not concerned with the individual pockets, and he 
further assumed that the variation in initial velocity 
required to predict the outcome would be immea- 
surable. It is while describing the game of roulette 
that Poincare introduces the concept of sensitivity 
to initial conditions, which is now a cornerstone of 
modern chaos theorjff^l. 



2 



16 2 



FIG. 1. The European roulette wheel. In the left panel one can see a portion of the rotating roulette wheel and 
surrounding fixed track. The ball has come to rest in the green pocket. Although the motion of the wheel and 
the ball (in the outer track) are simple and linear, one can see the addition of several metal deflectors on the stator 
(that is the fixed frame on which the rotating wheel sits). The sharp frets between pockets also introduce strong 
nonlinearity as the ball slows and bounces between pockets. The panel on the right depicts the arrangement of the 
number to 36 and the colouring red and black. 



A general procedure for predicting the outcome 
of a roulette spin, and an assessment of its utility 
was described by Edward Thorp in a 1969 publica- 
tion for the Review of the International Statistical 
Institute. In that paper. Thorp describes the two 
basic methods of prediction. He observes (as others 
have done later) that by minimising systematic bias 
in the wheel, the casinos achieve a mechanical per- 
fection that can then be exploited using determinis- 
tic prediction schemes. He describes two determin- 
istic prediction schemes (or rather two variants on 
the same scheme). If the roulette wheel is not per- 
fectly level (a tilt of 0.2° was apparently sufficient — 
we verified that this is indeed more than sufficient) 
then there effectively is a large region of the frame 
from which the ball will not fall onto the spinning 
wheel. By studying Las Vegas wheels he observes 
this condition is meet in approximately one third of 
wheels. He claims that in such cases it is possible 
to garner a expectation of -t-15%, which increased 
to -1-44% with the aid of a 'pocket-sized' computer. 
Some time later, Thorp revealed that his collabo- 
rator in this endeavour was Claude Shannorfi^, the 
founding father of information theorjffS]_ 

In his 1967 boolP the mathematician Richard A. 
Epstein describes his earlier (undated) experiments 



with a private roulette wheel. By measuring the an- 
gular velocity of the ball relative to the wheel he was 
able to predict correctly the half of the wheel into 
which the ball would fall. Importantly, he noted that 
the initial velocity (momentum) of the ball was not 
critical. Moreover, the problem is simply one of pre- 
dicting when the ball will leave the outer (fixed rim) 
as this will always occur at a fixed velocity. How- 
ever, a lack of sufficient computing resources meant 
that his experiments were not done in real time, and 
certainly not attempted within a casino. 

Subsequent to, and inspired by, the work of Thorp 
and Shannon, another widely described attempt to 
beat the casinos of Las Vegas was made in 1977-1978 
by Doyne Farmer, Norman Packard and colleague^ 
It is supposed that Thorp's 1969 paper had let the 
cat out of the bag regarding profitable betting on 
roulette. However, despite the assertions of Bas^, 
Thorp's papeil^is not mathematically detailed (there 
is in fact no equations given in the description of 
roulette). Thorp is sufficiently detailed to leave the 
reader in no doubt that the scheme could work, but 
also vague enough so that one could not replicate 
his effort without considerable knowledge and skill. 
Farmer, Packard and colleagues implemented the 
system on a 6502 microprocessor hidden in a shoe. 
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and proceeded to apply their method to the various 
casinos of the Las Vegas Strip. The exploits of this 
group are described in detail in Basa^. The same 
group of physicists went on to apply their skills to 
the study of chaotic dynamical system^i^l and also 
for profitable trading on the financial market^^SI Jn 
Farmer and Sidorowich's landmark paper on predict- 
ing chaotic time serie^^D the authors attribute the 
inspiration for that work to their earlier efforts to 
beat the game of roulette. 

Less exalted individuals have also been employ- 
ing similar schemes, in some cases fairly recently. 
In 2004, the BBC carried the report of three gam- 
blers (described only as "a Hungarian woman and 
two Serbian men"i2ll) arrested by police after win- 
ning £1, 300, 000 at the Ritz Casino in London. The 
trio had apparently been using a laser scanner and 
their mobile phones to predict the likely resting place 
of the ball. Happily, for the trio but not the casino, 
they were judged to have broken no laws and allowed 
to keep their winningJ^S. The scheme we describe in 
Section [H] and implement in Section |HI| is certainly 
compatible with the equipment and results reported 
in this case. In Section |IV] we conclude with some re- 
marks concerning the practicality of applying these 
methods in a modern casino, and what steps a casino 
could take (or perhaps have taken) to circumvent 
these exploits. A preliminary version of these results 
was presented at a conference in MacacF^. An inde- 
pendent and much more detailed model of dynamics 
of the roulette wheel is discussed in Strzalko et al^. 
Since our preliminary publicatiorl^, private commu- 
nication with several individuals indicates that these 
methods have now progressed to the point of at least 
four instances of independent in situ field trials. 





FIG. 2. The dynamic model of ball and wheel. On 

the left we show a top view of the roulette wheel (shaded 
region) and the stator (outer circles). The ball is moving 
on the stator with instantaneous position (r, 9) while the 
wheel is rotating with angular velocity tp (Note that the 
direction of the arrows here are for illustration only, the 
analysis in the text assume the same convention, clock- 
wise positive, for both ball and wheel). The deflectors 
on the stator are modelled as a circle, concentric with 
the wheel, of radius rdcfi. On the right we show a cross 
section and examination of the forces acting on the ball 
in the incline plane of the stator. The angle a is the 
incline of the stator, m is the mass of the ball, ac is the 
radial acceleration of the ball, and g is gravity. 



necessary. Moreover, it is messy. Finally, we suppose 
that the incline of the stator to the horizontal is a 
constant a. This situation, together with a balance 
of forces is depicted in figure [2j We will first consider 
the ideal case of a level table, and then in section II B 
show how this condition is in fact critical. 



II. A MODEL FOR ROULETTE 



We now describe our basic model of the motion 
of the roulette wheel and ball. Let (r, 6*) denote the 
position of the ball in polar co-ordinates, and let 
if denote the angular position of the wheel (say, the 
angular position of the centre of the green pocket) . 
We will model the ball as a single point and so let 
Ti-iin be the farthest radial position of that point (i.e. 
the radial position of the centre of the ball when the 
ball is spinning with high velocity in the rim of the 
wheel). Similarly, let rdefl be the radial distance to 
the location of the metal deflectors on the stator. For 
now, we will assume that = (that is, there 

are deflectors evenly distributed around the stator 
at constant radius rdcfi < t"). The extension to the 
more precise case is obvious, but, as we will see, not 



A. Level table 

For a given initial motion of ball (r, 9, 9, 9)t=o and 
wheel {if, if, (p)t=Q our aim is to determine the time 
idcfl at which r = rdcfi . After launch the motion of 
the ball will pass through two distinct states which 
we further divide into four cases: (i) with sufficient 
momentum it will remain in the rim, constrained by 
the fixed edge of the stator; (ii) at some point the 
momentum drops and the ball leaves the rim; (iii) 
the ball will gradually loose momentum while travel- 
ling on the stator as 9 drops, so will r; and (iv) even- 
tually r = rdefi at some time idefi- At time t = idcff 
we assume that the ball hits a deflector on the stator 
and drops onto the (still spinning) wheel. Of course, 
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the deflectors are discrete and located only at spe- 
cific points around the edge of the wheel. While it 
is possible, and fairly straightforward to incorporate 
the exact position (and more importantly, the ori- 
entation) of each deflector, we have not done this. 
Instead, we model the deflectors at a constant radial 
distance around the entire rim. The exact position 
of the wheel when the ball reaches the deflectors will 
be random but will depend only on (^(tdofl) — i-e. 
depending on where the actual deflectors are when 
the ball first comes within range, the radial distance 
until the ball actually deflects will be uniformly dis- 
tributed on the interval [0, 27r/iVdofl] where N^cR is 
the number of deflectors. 



(i) Ball rotates in the rim 

While traveling in the rim r is constant and the 
ball has angular velocity 9. Hence, the radial ac- 
celeration of the ball is Oc = ^ = lirO)'^ — rfP 
where v is the speed of the ball. During this pe- 
riod of motion, we suppose that r is constant and 
that Q decays only due to constant rolling friction: 
hence f = and Q = 0(0), a constant. This phase of 
motion will continue provided the centripetal force 
of the ball on the rim exceeds the force of gravity 
mac cos a > mg sin a (to is the mass of the ball) . 
Hence, at this stage 



■ „ g 

3 > - tana. 



(1) 



(ii) Ball leaves the rim 

Gradually the speed on the ball decays until even- 
tually 0^ = ^tanof. Given the initial acceleration 

0(0), velocity ^(0) and position 9{0), it is trivial to 
compute the time at which the ball leaves the rim, 
trim to be 
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^■(0) 



^(0) 



tan a 



(2) 



To do so, we assume that the angular acceleration 
is constant and so the angular velocity at any time 
is given by 9{t) = ^(0) + 9{Q)t and substitute into 
equation Q. That is, we are assuming that the 
force acting on the ball is independent of velocity — 
this is a simplifying assumption for the naive model 
we describe here, more sophisticated alternatives are 
possible, but in all cases this will involve the estima- 
tion of additional parameters. The position at which 



the ball leaves the rim is given by 

(ftana)-^(0)2 



271- 



e{o) . 

261(0) 

where | • |27r denotes modulo 27r. 
(iii) Ball rotates freely on the stator 



After leaving the rim the ball will continue (in 
practise, for only a short while) to rotate freely on 
the stator until it eventually reaches the various de- 
flectors at r = Tdcfl. The angular velocity continues 
to be governed by 



but 



that 



9{t) = 0(0) + 9{0% 



rO < gtana 



the radial position is going to gradually decrease too. 
The difference between the force of gravity mg sin a 
and the (lesser) centripetal force mr9^ cos a provides 
inward acceleration of the ball 



' cos a — gsma. 



(3) 



Integrating ([3| yields the position of the ball on the 
stator. 



(iv) Ball reaches the deflectors 

Finally, we flnd the time t = tdefl for which 
r{t), computed as the definite second integral of 
([3]), is equal to rdofl- We can then compute the in- 
stantaneous angular position of the ball 0(tdcfl) = 
61(0) + 0(O)tdcfi + ^9i0)tlcfi and the wheel (/^(tdefl) = 
<p(0) -|- ip{0)tdcR + |'^(0)i^j,fjto give the salient value 



7 = \0{tdc{i) - 'P(tdcfi)| 



27r 



(4) 



denoting the angular location on the wheel directly 
below the point at which the ball strikes a deflec- 
tor. Assuming the constant distribution of deflectors 
around the rim, some (still to be estimated) distribu- 
tion of resting place of the ball will depend only on 
that value 7. Note that, although we have described 
{9, 9, 9)t=Q and ((/?, (/S, (p)t=o separately, it is possible 
to adopt the rotating frame of reference of the wheel 
and treat — Lp ss a. single variable. The analysis is 
equivalent, estimating the required parameters may 
become simpler. 
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FIG. 3. The case of the crooked table. The blue 
curve denotes the stability criterion (|6|, while the red 
solid line is the (approximate) trajectory of the ball with 
8{ti) + 2iT = d{t2) indicating two successive times of com- 
plete revolutions. The point at which the ball leaves the 
rim will therefore be the first intersection of this stability 
criterion and the trajectory. This will necessarily be in 
the region to the left of the point at which the ball's tra- 
jectory is tangent to ([6|, and this is highlighted in the 
figure as a green solid. Typically a crooked table will 
only be slightly crooked and hence this region will be 
close to 9 = but biased toward the approaching ball. 
The width of that region depends on 0{ti) — 6(t2), which 
in turn can be determined from (Im. 



We note that for a level table, each spin of the ball 
alters only the time spent in the rim, the ball will 
leave the rim of the stator with exactly the same 
velocity 9 each time. The descent from this point 
to the deflectors will therefore be identical. There 
will, in fact, be some characteristic duration which 
could be easily computed for a given table. Doing 
this would circumvent the need to integrate pi). 



B. The crooked table 

Suppose, now that the table is not perfectly level. 
This is the situation discussed and exploited by 
Thorp^. Without loss of generality (it is only an 
affine change of co-ordinates for any other orienta- 
tion) suppose that the table is tilted by an angle S 
such that the origin </? = is the lowest point on the 
rim. Just as with the case of a level table, the time 



which the ball spends in the rim is variable and the 
time at which it leaves the rim depends on a stabil- 
ity criterion similar to ([ij. But now that the table 
is not level, that equilibrium becomes 

r9^ ^ gta.n {a + 5 cos 9). (5) 

If (5 = then it is clear that the distribution of 
angular positions for which this condition is first 
met will be uniform. Suppose instead that 6 > 0, 
then there is now a range of critical angular ve- 
locities G tan (a - 5), ^ tan (a -I- (5)]. Once 
9^ < ^ tan (a -I- 6) the position at which the ball 
leaves the rim will be dictated by the point of inter- 
section in {9, ^)-space of 



9= ^^tiin{a + Scos9) (6) 

and the ball trajectory as a function of t (modulo 
27r) 

9it) = 9(0) + 9{0)t (7) 

9it) = 9i0)+9{0)t+^9{0)t^. (8) 

If the angular velocity of the ball is large enough then 
the ball will leave the rim at some point on the half 
circle prior to the low point {ip = 0). Moreover, sup- 
pose that in one revolution (i.e. 9{ti) + 2TT = 9{t2) ), 
the velocity changes by 9{ti) — 9{t2)- Furthermore, 
suppose that this is the first revolution during which 
^2 < a tan (a -I- (5) (that is, 9{tif > f tan (a -I- 5) 
but 9{t2)'^ < Ha.n{a + S)). Then, the point at 
which the ball will leave the rim will (in {9, ^)-space) 
be the intersection of (|6| and 

9 ^ 9ih) - ^ [9{t2) - 9ih)) 9. (9) 

The situation is depicted in figure |3] One can ex- 
pect for a tilted roulette wheel, the ball will sys- 
tematically favour leaving the rim on one half of 
the wheel. Moreover, to a good approximation, the 
point at which the ball will leave the rim follows a 
uniform distribution over significantly less than half 
the wheel circumference. In this situation, the prob- 
lem of predicting the final resting place is signifi- 
cantly simplified to the problem of predicting the 
position of the wheel at the time the ball leaves the 
rim. 

We will pursue this particular case no further here. 
The situation ([5| may be considered as a generali- 
sation of the ideal 6 = case. This generalisation 
makes the task of prediction significantly easier, but 
we will continue to work under the assumption that 
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the casino will be doing its utmost to avoid the prob- 
lems of an improperly levelled wheel. Moreover, this 
generalisation is messy, but otherwise uninteresting. 
In the next section we consider the problem of im- 
plementing a prediction scheme for a perfectly level 
wheel. 



III. EXPERIMENTAL RESULTS 

In Sec. HI] we introduces the basic mathemati- 
cal model which we will utilise for the prediction of 
the trajectory of the ball within the roulette wheel. 
We ignore (or rather treat as essentially stochastic) 
the trajectory of the ball after hitting the deflectors 
— charting the distribution of final outcome from 
deflector to individual pocket in the roulette wheel 
is a tractable probabilistic problem, and one which 
for which we will sketch a solution later. However, 
the details are perhaps only of interest to the profes- 
sional gambler and not to most physicists. Hence, we 
are reduced to predicting the location of the wheel 
and the ball when the ball first reaches one of the de- 
flectors. The model described in Sec.,|ll]is sufRcient 
to achieve this — provided one has adequate mea- 
surements of the physical dimensions of the wheel 
and all initial positions, velocities and accelerations 
(as a further approximation we assume deceleration 
of both the ball and wheel to be constant over the 
interval which we predict. 

Hence, the problem of prediction is essentially 
two-fold. First, the various velocities must be es- 
timated accurately. Given these estimates it is a 
trivial problem to then determine the point at which 
the ball will intersect with one of the deflectors on 
the stator. Second, one must then have an estimate 
of the scatter imposed on the ball by both the deflec- 
tors and possible collision with the individual frets. 
To apply this method in situ, one has the further 
complication of estimating the parameters r, rdcfl. 
Trim, a and possibly S without attracting undue at- 
tention. We will ignore this additional complication 
as it is essentially a problem of data collection and 
statistical estimation. Rather, we will assume that 
these quantities can be reliably estimated and re- 
strict our attention to the problem of prediction of 
the motion. To estimate the relevant positions, ve- 
locities and accelerations {9, 6, 6, if, ip, </3)t=o (or per- 
haps just [9 — if,6 — Lp,9~ V3)t=o) we employ two 
distinct techniques. 

In the following subsections we describe these 
methods. In Sec. IIII Al we introduce a manual 



measurement scheme, and in Sec. IIIB we describe 



our implementation of a more sophisticated digital 
system. The purpose of Sec. HI A is to demon- 
strate that a rather simple "clicker" type of device — 
along the lines of that utilized by the Doyne Farmer, 
Norman Packard and coUaboratorJi] — can be em- 
ployed to make sufficiently accurate measurements. 
Nonetheless, this system is far from optimal: we 
conduct only limited experiments with this appa- 
ratus: sufRcient to demonstrate that, in principle, 
the method could work. In Sec. IIII Bl we describe 
a more sophisticated system. This system relies on 
a digital camera mounted directly above a roulette 
wheel and is therefore unlikely to be employed in 
practice (although alternative, more subtle, devices 
could be imagined). Nonetheless, our aim here is to 
demonstrate how well this system could work in an 
optimal environment. 

Of course, the degree to which the model in Sec. 
[IT] is able to provide a useful prediction will depend 
critically on how well the parameters are estimated. 
Sensitivity analysis shows that predicted outcome 
(Eq. [4]) depends only linearly or quadratically (in 
the case of physical parameters of the wheel) on 
our initial estimates. More important however, and 
more difficult to estimate, is to what extent each of 
these parameters can be reliably estimated. For this 
reason we take a strictly experimental approach and 
show that even with the various imperfections inher- 
ent in experimental measurement, and in our model, 
sufficiently accurate predictions are realizable. 



A. A manual implementation 

Our first approach is to simply record the time 
at which ball and wheel pass a fixed point. This is 
a simple approach (probably that used in the early 
attempts to beat the wheels of Las Vegas) and is 
trivial to implement on a laptop computer, personal 
digital assistant, embedded system, or even a mo- 
bile phon026ii27i28j^ Q^j, results, depicted in figure 
|4] illustrate that the measurements, although noisy, 
are feasible. The noise introduced in this measure- 
ment is probably largely due to the lack of physical 
hand-eye co-ordination of the first author. Figure |4] 
serves only to demonstrate that, from measurements 
of successive revolutions T{i) and T{i + 1) the rela- 
tionship between T{i) and T(i + 1) can be predicted 
with a fairly high degree of certainty (over several 
trials and with several different initial conditions. 
As expected, the dependence of T{i -I- 1) on T{i) is 
sub-linear. Hence, derivation of velocity and acceler- 
ation from these measurements should be relatively 
straightforward. Simple experiments with this con- 
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FIG. 4. Hand-measurement of ball and wheel ve- 
locity for prediction. From two spins of the wheel, 
and 20 successive spins of the ball we logged the time 
(in seconds) T{i) for successive passes past a given point 
(T(i) against T{i + 1)). The measurements T(i) and 
T{i + 1) are the timings of successive revolutions — di- 
rect measurements of the angular velocity observed over 
one complete rotation. To provide the simplest and most 
direct indication that handheld measurements of this 
quantity are accurate, we indicate in this figure a de- 
terministic relationship between these quantities. From 
this relationship one can determine the angular deceler- 
ation. The red points depict these times for the wheel, 
the blue points are for the ball. A single trial of both ball 
and wheel is randomly highlighted with crosses (super- 
imposed). The inset is an enlargement of the detail in 
the lower left corner. Both the noise and the determin- 
ism of this method are evident. In particular, the wheel 
velocity is relatively easy to calculate and decays slowly, 
in contrast the ball decays faster and is more difficult to 
measure. 

figuration indicate that it is possible to accurately 
predict the correct half of the wheel in which the 
ball will come to rest (statistical results are given in 
the caption of Fig. |4]) . 

Using these (admittedly noisy) measurements we 
were able to successfully predict the half of the wheel 
in which the ball would stop in 13 of 22 trials (ran- 
dom with p < 0.15), yielding an expected return of 
36/18 X 13/22 - 1 = +18%. This trial run included 
predicting the precise location in which the ball 
landed on three occasions (random with p < 0.02). 
Quoted p-values are computed against the null hy- 
pothesis of a binomial distribution of n trials with 
probability of success determine by the fraction of 
the total circumference corresponding to the target 
range — i.e. the probability of landing by chance in 



one of the target pockets. 

B. Automated digital image capture 

Alternatively, we employ a digital camera 
mounted directly above the wheel to accurately and 
instantaneously measure the various physical param- 
eters. This second approach is obviously a little 
more difficult to implement incognito. Here, we are 
more interested in determining how much of an edge 
can be achieved under ideal conditions, rather than 
the various implementation issues associated with 
realising this scheme for personal gain. In all our tri- 
als we use a regulation casino-grade roulette wheel (a 
32" "President Revolution" roulette wheel manufac- 
tured by Matsui Gaming Machine Co. Ltd., Tokyo). 
The wheel has 37 numbered slots (1 to 36 and 0) in 
the configuration shown in figure [T] and has a radius 
of 820 mm (spindle to rim). For the purposes of 
data collection we employ a Prosilica EC650C lEEE- 
1394 digital camera (1/3" CCD, 659x493 pixels at 
90 frames per second). Data collection software was 
written and coded in C++ using the OpenCV li- 
brary. 

The camera provides approximately (slightly less 
due to issues with data transfer) 90 images per sec- 
ond of the position of the roulette wheel and the 
ball. Artifacts in the image due to lighting had to 
be managed and filtered. From the resultant im- 
age the position of the wheel was easily determined 
by locating the only green pocket ("0") in the wheel, 
and the position of the ball was located by differenc- 
ing successive frames (searching for the ball shape or 
color was not sufficient due to the reflective surface 
of the wheel and ambient lighting conditions). 

From these time series of Cartesian coordinates for 
the position of both the wheel (green "0" pocket) 
and ball we computed the centre of rotation and 
hence derived angular position time series. Polyno- 
mial fits to these angular position data (modulo 27r) 
provided estimates of angular velocity and accelera- 
tion (deceleration). From this data we found that, 
for out apparatus, the acceleration terms where very 
close to being constant over the observation time pe- 
riod — and hence modeling the forces acting on the 
ball as constant provided a reasonable approxima- 
tion. With these parameters we directly applied the 
model of Section [IT] to predict the point at which the 
ball came into contact with the deflectors. 

Figure [5] illustrates the results from 700 trials of 
the prediction algorithm on independent rolls of a 
fair and level roulette wheel. The scatter plot of Fig. 
[5] provides only a crude estimation of variance over 
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FIG. 5. Predicting roulette. The plot depicts the re- 
sults of 700 trials of our automated image recognition 
software used to predict the outcome of independent 
spins of a roulette wheel. What we plot here is a his- 
togram in polar coordinates of the difference between the 
predicted and the actual outcome (the "Target" location, 
at the 12 o'clock position in this diagram, indicating that 
the prediction was correct). The length of each of the 
37 black bars denote the frequency with which predicted 
and actual outcome differed by exactly the correspond- 
ing angle. Dot, dot-dashed and solid (red) lines depict 
the corresponding 99.9%, 99% and 90% confidence in- 
tervals using the corresponding two-tailed binomial dis- 
tribution. Motion forward (i.e. ball continues to move 
in the same direction) is clockwise, motion backwards in 
anti-clockwise. From the 37 possible results there are 2 
instances outside the 99% confidence interval. There are 
7 instances outside the 90% confidence interval. 



the entire region of the wheel for a given prediction. 
A determined gambler could certainly extend this 
analysis with a more substantial data set relating 
to their particular wheel of interest. We only aim 
to show that certain non-random characteristics in 
the distribution of resting place will emerge and that 
these can then be used to further refine prediction. 

Nonetheless, several things are clear from Fig. [sj 
First, for most of the wheel, the probability of the 
ball landing in a particular pocket — relative to the 
predicted destination — does not differ significantly 
from chance: observed populations in 30 of 37 pock- 
ets is within the 90% confidence interval for a ran- 
dom process. Two particular pockets — the target 



pocket itself and a pocket approximately one-quarter 
of the wheel prior to the target pocket — occur with 
frequencies higher than and less than (respectively) 
the expected by chance: outside the 99% confidence 
interval. Hence, the predicted target pocket is a 
good indicator of eventual outcome and those pock- 
ets immediate prior to the target pocket (which the 
ball would need to bounce backwards to reach) are 
less likely. Finally, and rather speculatively, there is 
a relatively higher chance (although marginally sig- 
nificant) of the ball landing in one of the subsequent 
pockets — hence, suggesting that the best strategy 
may be to bet on the section of the wheel following 
the actual predicted destination. 

IV. EXPLOITS AND COUNTER-MEASURES 

The essence of the method presented here is to 
predict the location of the ball and wheel at the 
point when the ball will first come into contact with 
the deflectors. Hence, we only require knowledge 
of initial conditions of each aspect of the system 
(or more concisely, their relative positions, veloci- 
ties and accelerations). In addition to this, certain 
parameters derived from the physical dimensions of 
the wheel are required — these could either be es- 
timated directly, or inferred from observational tra- 
jectory data. Finally, we note that while anecdotal 
evidence suggests that (the height of the) frets plays 
an important role in the final resting place of the 
ball, this does not enter into our model of the more 
deterministic phase of the system dynamics. It will 
affect the distribution of final resting places — and 
hence this is going to depend rather sensitively on a 
particular wheel. 

We would like to draw two simple conclusions from 
this work. First, deterministic predictions of the 
outcome of a game of roulette can be made, and 
can probably be done in situ. Hence, the tales of 
various exploits in this arena are likely to be based 
on fact. Second, the margin for profit is quite slim. 
Minor manipulation with the frictional resistance 
or level of the wheel and/or the manner in which 
the croupier actually plays the ball (the force with 
which the ball is rolled and the effect, for exam- 
ple, of axial spin of the ball) have not been explored 
here and would likely affect the results significantly. 
Hence, for the casino the news is mostly good — 
minor adjustments will ameliorate the advantage of 
the physicist-gambler. For the gambler, one can rest 
assured that the game is on some level predictable 
and therefore inherently honest. 

Of course, the model we have used here is ex- 
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tremely simple. In Strzalko et alW^ much more so- 
phisticated modehng methodologies have been in- 
dependently developed and presented. Certainly, 
since the entire system is a physical dynamical sys- 
tem, computational modeling of the entire system 
may provide an even greater advantagd^. Nonethe- 
less, the methods presented in this paper would cer- 
tainly be within the capabilities of a 1970's "shoe- 
computer" . 
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